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Abstract 

There is a need to go beyond the narrow resonance approximation for QCD sum-rule 
channels which are likely to exhibit sensitivity to broad resonance structures. We first discuss 
how the first two Laplace sum rules are altered when one goes beyond the narrow resonance 
approximation to include possible subcontinuum resonances with nonzero widths. We then show 
that the corresponding first two finite energy sum rules are insensitive to the widths of such 
resonances, provided their peaks are symmetric and entirely below the continuum threshold. We 
also discuss the reduced sensitivity of the first two finite energy sum rules to higher dimensional 
condensates, and show these sum rules to be insensitive to dimension > 6 condensates containing 
at least one qq pair. We extract the direct single-instanton contribution to the Fj sum rule for the 
longitudinal component of the axial-vector correlation function from the known single-instanton 
contribution to the lowest Laplace sum rule for the pseudoscalar channel. Finally, we demonstrate 
how inclusion of this instanton contribution to the finite-energy sum rule leads to both a lighter 
quark mass and to more phenomenologically reasonable higher-mass-resonance contributions 
within the pseudoscalar channel. 



I. Introduction: Nonzero Resonance Widths and QCD Laplace Sum-Rules 

Laplace Sum-Rule Methodology in the Narrow Resonance Approximation 

Hadron properties can be extracted by relating phenomenological and field-theoretical 
expressions for integrals over appropriately chosen current-correlation functions, integrals which 
we denote as QCD sum rules [1]. The phenomenological expressions are generally extracted 
via the narrow resonance approximation. In the narrow resonance approximation, hadronic 
contributions to the imaginary part of current-current correlation functions are proportional to 5- 
functions at the resonance mass, 

Im[n\s)] = E 7ig,5(s - m^) + 0(s - So) Im[nP(s)], (1) 

The summation in (1) is over all resonances r in the channel under consideration [i.e., whose 
quantum numbers are consistent with the choice of currents in the current correlation function] 
such that mj is less than Sq. Above this hadron-continuum threshold, the hadronic contribution 
n''(s) to the correlation function is assumed to be the same as the contribution nP(s) from 
perturbative QCD, as is evident from (1). 

The hadronic sub-continuum (h) contribution to the k* Laplace sum rule, corresponding 
to the transform of the appropriate portions of (1), is defined to be 

oo 

R^(T) = / ds (1/7U) Im[n^s) - nP(s) 0(s - Sq)] e"'^ (2) 


In the narrow resonance approximation (F — > 0), we see from (1) that 
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Lim R^(T) = E mf expL-m^T], (3) 

r^o 

an expression in which contributions from more-massive resonances are exponentially suppressed. 
Note from (3) that R\(x) > m^ Ro(x), where m^ denotes the mass of the lowest-lying resonance 
in the channel. Consequently, R'i'(x)/RQ(x) is bounded from below by m^. Standard QCD sum- 
rule methodology involves minimizing this ratio [or its field-theoretical analogue] with respect 
to T in order to determine a value of m^ [2]. The sum rule R^Ct) corresponds to the following 
field-theoretical contribution from perturbative-QCD and nonperturbative (np) QCD-vacuum 
effects: 

So 

rQCD(^) ^ / (yjj) Im[nP(s)] s'' e-'^ 



+ (-a/aT)'^{(i/T) sf-i[-dn"p(s)/dQ2]}. (4) 

In equation (4), = -s, and n"P(s) represents all correlation-function contributions from QCD- 
vacuum condensates as well as additional finite-correlation length contributions from the 
instanton background, i.e. "direct instanton contributions." The inverse Laplace transform in (4), 
corresponding to the Laplace-transform definition 

oo 

^Q.[f(T)] = / dxf(x)e-Q^^ (5) 



is utilized to take advantage of the operator-product expansion of 11"^ in inverse powers of Q ^ , 
and is easily understood via dispersion-relation methodology. Noting first that the singularities 
of n'*(s) [hadron poles and kinematic production-threshold branch cuts] must lie on the positive 
real s-axis, one finds that 



(1/x) se-^[-dn"P(s)/dQ2] 



(1/x) S£^^[(1/7C)/ ds Im[n"P(s)]/(s+Q2)2] 

(6) 
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As is evident from (5), 9E^^[l/(s+Q2)2] = xe which, upon substitution into (6) and (4), leads 
to a result consistent with duality between QCD [nP(s) + n"P(s)] and phenomenological hadronic 
physics [n\s)]: 



R^c°(T) + / ds (1/71) Im[nP(s)] s'' e 



-St 



/ ds (l/7t) Im[nP(s) + n"P(s)] Q-'^ 

(7) 



Duality between R^^'^(t) and then follows via comparison of equations (7) and (2). We 
then find that the lowest lying resonance can be determined via the relationship: 

Min[R?'=^°(T)/R^^^(T)] > m] (8) 



over an appropriate range of T [Sq > x » Aq^jj] 



Laplace Sum-Rule Width Corrections to the Lowest-Lying Resonance Mass 

There is a need to go beyond the narrow resonance approximation if QCD sum rules 
exhibit sensitivity to resonance structures with non-zero widths. Such structures can not always 
be absorbed in the sum-rule continuum—even the lowest hadronic resonances may have 
substantial widths. For example, theoretical arguments exist [3,4] for the first pion-excitation 
to have a mass below 1 GeV, a floor for any reasonable estimate of the continuum threshold 
above which perturbative and hadronic QCD should coincide. Even if the first pion excitation 
state is identified with the 11(1300) resonance, whose mass pole is still likely to be below the 
continuum threshold, the width of this resonance may be as large as 600 MeV [5]. Moreover, 
the lowest isoscalar 0"^ meson, if it exists at all [6], may have a width even larger than its 
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mass [7], though arguments for a a(550)-resonance with a somewhat more moderate width 
have also been recently advanced [8]. 



To gain qualitative insight into how nonzero resonance widths can effect QCD sum rule 
calculations, we can replace the 5-function of a resonance contribution to (1) with a rectangular 
pulse of unit area: 



Equation (9) defines a rectangular pulse centred at s = m^ with full- width As = 2mr and height 
l/(2mr). 

Let us consider how such an approximation to a lowest-lying resonance alters a QCD 
Laplace sum-rule determination of that resonance's mass. We assume in the spirit of the original 
formulation of QCD sum rules [1] that all but the lowest-lying (C) resonance is absorbed in the 
continuum. If we replace the delta function for the lowest-lying resonance with the pulse 

(?jn(s,r), we find from (2) that 



5(s-m2) ^J.s,r) = [0(s-m2+mr) - 0(s-m2-mr)]/2mr. 



(9) 



g, /°dsPJs,ne-^^ = 



Ao(m,r,T), 



(10) 







f ds P^(s,r) s e-^^ = 



gfiin^ e 



Ai(m,r,T), 



(11) 







with the functions Aq j found from explicit evaluation of the integrals in (10) and (11): 



AQ(m,r,x) = sinh(mrx)/(mrx) 



(12) 



Ai(m,r,x) = Ao(m,r,x)[l + l/(m2x)] - cosh(mrx)/(m2x). 



(13) 
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The results (10-13) assume that P < m, Sq ^ + mr, so that the integration includes the entire 
resonance peak. Note also that 

Lim A(ji(m,r,T) = 1, 
r^O 

consistent with the 5-function limit of a square pulse of infinitesimal width. 

We see immediately from (10) and (11) that 

^ R?(x) Ao(m,r,x ) 

rS(t) Ai(m,r,T) (14) 

Since rJ(t)/Ro(t) corresponds to R'^^^/R^*"^^ by duality, and since this latter ratio corresponds 
to m^ in the narrow resonance approximation (r=0), one can show from (12) and (13) that finite 
width effects will increase the masses of lowest-lying resonances extracted via Laplace sum 
rules: 

= [m2]r=o " [Ao(m,r,T)/Ai(m,r,T)] 
= [m2]r^o-[l+r^x/3 + 0(r4)]. (15) 

This result should properly be regarded as a lower-bound on the magnitude of width 
contributions to R/Rq- The width F appearing in (15) is, at present, defined via the rectangular 
pulse (9); it cannot be understood to correspond to a Breit-Wigner resonance width. Indeed, the 
narrow resonance approximation follows from the narrow-width limit of the Breit-Wigner 
resonance shape: 

Im[n^(s)] = Im[-g/(s-m2+imn] 

= g^r/[(s-m2)2+ m^r^] 7tg^5(s-m2). (16) 

r^o 

If a resonance has a Breit-Wigner shape, then half of the total area of the resonance-peak 
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(considered as a function of s) is included in the range - mP < s < + mP. Since the unit- 
area rectangular pulse (9) has all of its area included in the range m^- mP < s < m2+ mP, the 
result (15) is based upon a narrower pulse than the Breit-Wigner pulse of equivalent P, and is 
therefore likely to be an underestimate of the contribution of a Breit-Wigner resonance-width P 
to the Laplace sum-rule determination of m^. 

For Laplace sum rules, a more quantitative estimate of resonance-width effects could be 
obtained by replacing the delta-functions in (1) with the Breit-Wigner peaks (16), and then 
substituting into the Laplace sum-rule definition (2). However, such an approach is subject to 
ambiguity. The Breit-Wigner shape has an infinite tail, and significant portions of that tail may 
extend above the continuum threshold Sq or below the s=0 boundary into Euclidean momenta. 
Truncating such contributions would artificially exclude some of the integrated resonance peak. 
On the other hand, including such contributions leads to methodological contradictions with 
hadron physics. The Breit-Wigner shape, which itself stems from a linear approximation, can be 
modified for broad widths so as to vanish at s = [9]. Even with such a modification, post- 
continuum contributions from the Breit-Wigner tail, whether included or truncated away, can be 
genuinely substantial for resonances with widths in excess of 100 MeV, and can be a source of 
theoretical uncertainty in Laplace sum-rule analyses of broad sub-continuum resonances. 

Such uncertainty may be understood as a limitation on Laplace sum-rule methodology 
itself, particularly for channels in which more than one resonance lies below the continuum 
threshold. Non-lowest-lying resonances are expected to be less stable, and consequently, to be 

substantially broader than lowest- lying resonances. The 1=1 pseudoscalar channel has already 
been mentioned as an example of such a channel, and is discussed in the final two sections of 
this paper. ^ 



A QCD sum-rule treatment of the 1 = scalar channel is even more likely to be 
problematical, in that this channel may be sensitive to not only the controversial sigma 
[fQ(400-1200)], but also as many as three other subcontinuum resonances [fQ(980), 
fQ(1370), and fQ(1500)], with fQ(1370) being a broad resonance. 
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In the section that follows, we first discuss how finite-energy sum rules can alleviate the 
resonance-width ambiguities described above. Unlike the case for Laplace sum rules, we 
demonstrate that the contribution of a non-narrow resonance to the first two finite-energy sum 
rules (Fq and Fj) is independent of that resonance's width, provided the resonance peak is both 
symmetric and entirely below the continuum threshold. 

These same two sum rules are also less sensitive to higher dimensional condensates than 
corresponding Laplace sum rules, thereby lowering the number of nonperturbative QCD-vacuum 
parameters required to enter a sum-rule analysis. In Section III, we demonstrate that Fq and F^ 
are essentially decoupled (to leading order in a^) from dimension > 6 QCD-vacuum condensates 
containing one or more qq pairs. For the specific channel pertinent to pseudoscalar resonances, 
this suppression is shown to occur even for dimension > 4; the dimension- 6 condensate (a5(qq) ^ ) 
is seen (by virtue of a third order pole at =0) not to enter Fq j, even though it is known to 
enter the corresponding Laplace sum rules Rq Insensitivity of the first two finite-energy sum 
rules to multiple gluon condensates [condensates which do not have any qq pairs] of dimension 
> 6 is discussed at length in Section IV, including the operator-mixing that serves to suppress 
the dimension-6 gluon condensate. 

We then focus on QCD sum rules for the I = 1 pseudoscalar mesons, as noted earlier, as 
an example of a "problem" channel with broad subcontinuum resonance contributions. This 
channel has long been understood to be subject to instanton contributions. In Section V, we 
extract the direct single-instanton contribution to the finite-energy sum rule Fj for this channel. 
In Section VI, we relate F^ to the phenomenology of this channel, which has been used by others 
[10,11,12] to obtain a fairly large lower bound on the light quark mass (m^ + m^ > 20 ± 
5 MeV). We demonstrate how direct single-instanton contributions to Fj not only lower the 
bound for the estimated quark mass, but also lower the estimated contribution of higher-mass 
resonances to the sum rule to values of r^ [= (FiMj/fj^m^) ^ ] more consistent with present 
phenomenological expectations. 
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Our manuscript also possesses four detailed appendices pertinent to a full methodological 
understanding of finite energy sum rules. In Appendix A an exact expression for the imaginary 
part of a function that arises in closed fermion loop contributions to correlation functions [i.e. in 
purely perturbative and multiple-gluon condensate contributions] is extracted. In Appendices B 
and C, the two-gluon condensate contributions to Fq and Fj are respectively calculated from that 
condensate's exact (as opposed to leading-order in m ) one-loop contribution to the longitudinal 
component of the axial vector correlation function. This contribution is shown to arise entirely 
from a net branch singularity for s > 4m 2. The absence of net pole contributions at s = 0, as 
well as the cancellation of infrared singularities arising from integration of the exact expression 
along the branch cut against those arising from integration around the branch-cut terminus at s 
= 4m 2 , are also demonstrated explicitly. All of these results (including the singularity structure 
described above) are applicable to the two-gluon condensate contributions to the finite energy 
sum rules associated with the scalar, vector and the transverse component of the axial vector 
correlation functions. These sum rule contributions are itemized in Appendix D. The explicit 
cancellation of quark-mass singularities via operator mixing is also demonstrated for channels in 
which they naively occur. 

II. Nonzero Resonance Widths and Finite-Energy Sum-Rules 

Finite Energy Sum Rules and Higher-Mass Resonances 

For a given current-correlation function n(s), the finite-energy sum-rules (FESR's) F,^(Sq) 
are defined here to be the integrals [13] 

F^(Sq) = (l/27ti) / dss''n(s), 

C(So) (17) 

where the contour C(so) is an open circle of radius Sq in the complex s-plane that does not cross 
the real s-axis [Fig. la]. The parameter Sq is understood to be the continuum threshold discussed 
in the previous section. For the hadronic contribution to the FESR's F^(Sq), the contour C(So) can 
be distorted into a line running below and above the physical singularities on the positive real 
s-axis [Fig. lb]: 
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F^(so) = (l/7t) / ds s'' Im[n\s)]. (18) 




In the narrow resonance approximation (1), one finds that 

F^(so) = E g,mf = E [F^(so)], (19) 

r r 

an expression that differs from (3) only in that higher-mass sub-continuum resonances are no 
longer exponentially suppressed. This is a positive feature of the FESR approach, if one is 
seeking to use sum rules to obtain information about such resonances. 

Insensitivity of Fq j to Symmetric-Peak Resonance Widths 

To examine finite width effects, let us first replace the delta-functions of (1) with the 
finite-width rectangular pulses (9). As long as Sq > m^ + 2m^p the contribution of such a pulse 
to Fq is clearly the same as that of a delta-function, since Fq is sensitive only to peak-area: 

[FS(so)], ^ /"dsg,P^(s,r,) = g, (20) 



Remarkably, the Fj sum-rule is also insensitive to the width of the rectangular pulse: 



[F\(sq)1 ^ / ds s g,P^(s,r,) = [g/(2m^,)]r ds s = g^? (21) 

J 2 



The final result of (21) is identical to the contribution to [F^(Sq)]j obtained from the narrow 
resonance approximation, with Pj^^(s,rj.) replaced by 5(s-m^). The results (20) and (21) are to be 
contrasted with the width-dependence exhibited in (10) and (11) for corresponding Laplace sum 
rules. 
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Moreover, the width-independence of the first two FESR's obtained above is not an 
artifact of the rectangular pulse approximation for non-zero width resonances. Any symmetric 
resonance peak ]R^(s) centred at m^ can be represented as a sum over variable- width unit-area 
rectangular pulses Pjjj(s,r') centred at s = m^ : 

r 

max 

E„(s) = / dr f(r) P^(s,r) (22) 




In (22), f(r') is just the weighting assigned to the unit-area rectangular pulse with width F'. For 
example, the Breit-Wigner shape (16) can be expressed in the form of (22) by converting a 
Riemann sum of infinitesimally thin pulses into an integral [Fig. 2]: 



mF/[(s-m2)2+m2F2] 

Lim {(2/n) E V(n/j)-l Pjs, V(n/j)-l F)} 

n— >oo j=l 



2 / dy V(l-y)/y 9Js, V(l-y)/y F) 





= 4Ff ^ll^^ P>,r) 



(23) 



Assuming the peak ^^(s) has an area normalized to %, consistent with Mm(s) — > % 5(s-m2) in the 
narrow resonance limit [e.g. eq. (16)], one finds that 



^0 '"max 

n= / ds E^(s) = / dr f(r) (24) 
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provided Sg > + mTj^^^.^ Consequently, one can use (21) and (24) to demonstrate that 
replacing factors of % d(s-m^) in (1) with M^(s) will not alter narrow-resonance approximation 
predictions (19) for Fg and F^: 



[FS(So)]r 



»0 

^ / ds (1/7C) g, [RJs)], = g. 



(25) 







[F5'(So)]r 



^ / ds S (l/7C)g,[]R^(s)], 








(26) 







Thus we see that the finite-energy sum-rules Fq and Fj are impervious to resonance-width 
effects, provided the resonance in question is a symmetric peak that is entirely below the 
continuum threshold Sq. Consequently, we observe that these sum rules are particularly well- 
suited for an analysis of non-lowest-lying subcontinuum resonances. As remarked earlier, the 
contributions (19) of such resonances to Fq and F^ are not exponentially suppressed, as is the case 
with Laplace sum-rules (3). Since such resonances are unstable, and therefore broad, it is 
phenomenologically useful that their contributions to Fq and Fj are unaffected by their decay 
widths. It should be noted this property does not apply to higher finite energy sum rules; sum 
rules Fy. do exhibit width dependence if k > 2. Consequently, we will restrict our discussion 
henceforth to the properties of Fq j sum rules. 



^ For the theoretical Breit-Wigner shape, problems associated with the infinite tail of the 
resonance extending past the continuum threshold Sq have already been mentioned. One 
approximation which truncates the Breit-Wigner shape symmetrically is to choose T^^^ 
= 2m. This approximation makes sense, however, only if F « 2m < Sq. 
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III. FESR Suppression of Higher-Dimensional Condensates with One or More qq-Pairs 

The General Case 

The operator-product expansion (OPE) for a dimension-2 two-current correlation function 
n(s) can be expressed at Euclidean momenta Q ^ = -s > in terms of QCD-vacuum condensates 
as follows: 

n(-Q2) = Cp(Q2) + Cqq(Q2)<mqqq> + CoaCQ^Xa^G^) + C^CQ^XqC-aq) 

+ CG3(Q2)<a,G3> + C(,<,.(Q2)<a,(qq)2> + ... (27) 
To leading order in a^, the OPE coefficients C„(Q of an n-dimensional condensate (Oj are of 
the general form 

C„(Q2) = E[Aj + Bjln(QV^l2)]mJ/Qi+"-2 (28) 

j 

To avoid mass singularities, the index j is restricted to zero and even positive integers if n is 
even, and to odd positive integers if n is odd. To leading order in a^, contributions to (27) from 
condensates containing at least one fermion-antifermion pair necessarily correspond to diagrams 
with broken loops [Fig. 3], and for such diagrams Bp = 0; logarithms from integrations over 
closed-loop momenta do not occur. For example, to leading order in a^, the (m^qq) contribution 
[Fig. 3a] to the longitudinal component n^(s) of the axial-vector current correlation function, 

(g^^- p^yp2)nT(p2) + (p^yp2)nHp^) ^ i/d^x e^P-" <0|T j^5(x) j,5(0)|0>, (29) 
IJn5 = uy^Ysd] is given by [14] 

C-;(Q=^) = (2/m2)[l - (1 + 4m2/Q2)^] = -4/Q2 + Am^lQ' - Sm^/Q^ + ... (30) 

If in (28) Bj = for all j, the definition (17) implies that the FESR's Fq and Fj are (respectively) 
sensitive only to first and second order poles at =0: 
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[Fo(So)]qq= (l/2ra)/ ds Cq^C-sXmqqq) = -4<mqqq> (31) 

C(So) 

[Fi(So)]qq= (l/2ra)/ ds s C-L(-sXmqqq> = -4m2 (m^qq) (32) 

C(So) 

Thus, if Cjj(Q 2 ), the OPE coefficient of a condensate (0„>, is restricted to inverse powers of Q ^ , 
then n must be less than or equal to 6 for that condensate to contribute to Fq or Fj. If n > 6, 
then n+j-2 > 6 and the leading OPE contribution to (28) is at least a third order pole at = 0, 
which cannot contribute to Fq or F^ 

Higher-Dimensional Fermionic Condensates in the Pseudoscalar Channel 

For the particular case of the longitudinal component (L) of the axial-vector correlation 
function, which is coupled to pion-resonance states, there is an additional chiral symmetry 
constraint that C^(Q ^ ) — > as m^ ^ 0, in which case j > 1 for all coefficients of condensates 
that fail to vanish in the chiral limit. As a consequence, one can show to leading order in that 
the n = 6 condensate {QLj^qc{)^) cannot contribute to Fq or Fj, as its leading contribution is 
necessarily a third-order pole at = [1,10]: 

C(^,)^(Q ^ ) = -4487U m\ aJ2iqf> + 0(mJ/Q8). (33) 

Similarly, Fq and Fj are found to be insensitive to the (n = 5) mixed condensate (qG-oq). The 
relevant contribution to the longitudinal component of the axial-vector correlator is also seen to 
involve only third-and-higher order poles at = [15]: 

Cii(Q') = -(l-v)'/2m3v = 4m3/Q6 - 20m5/Q8 +..., (34) 



V = (1 + 4m^/Q2)^ (35) 
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Thus the leading contributions to Fq and sum rules in this channel do not involve any 
condensates with quark-antiquark pairs except (m^qq). The Fq and Fj sum rules in other channels 
can also involve the n=5 mixed condensate (qG-aq> and the n=6 condensate {a^(qq)^) [we are 
assuming vacuum-saturation], but no other condensates containing quark-antiquark pairs, as all 
other such condensates are of dimension greater than 6. 



IV. Purely Gluonic Contributions to and Fj 

Purely Perturbative Gluon-Loop Contributions 

For two-current correlation functions, the suppression of leading-order contributions from 
n > 6 condensates applies only to those operators whose leading contribution in does not 
involve a closed perturbative loop. However, all condensates involving gluons necessarily are 
generated from the closed-loop vacuum polarization diagram [Fig. 4], and such diagrams are 
characterized by nonzero coefficients Bj in the OPE expansion (28). The contribution of such 
logarithmic terms in (28) to the FESRs Fq and F^ can be obtained from the general relation [Q^ 
= -s; D is an integer] 

j ds ln(Q2)/(Q2)D = -2i7t(-l)° s^-°/(l-D); D>2 (36) 
C(So) 

However, a more precise evaluation of the contributions of closed-loop OPE coefficients 
necessarily involves the one-loop momentum integral^ 

1 

X(v) = (l/v2) / dx ln[l - s x(l-x)/m^ - ilel] +2/v\ (37) 



V = (1 - 4m^/s)''. (38) 

SI 



^ In eqs. (37-9) we are utilizing the notation of ref. [15]. 
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For Euclidean momenta (s < 0), one finds X(v) to be the real function 

X(v) = (1/v) ln[(l+v)/(v-l)]. (39) 

For Minkowskian momenta (s > 0), X(v) develops an imaginary part above the quark- antiquark 
kinematic production threshold, as discussed in Appendix A: 

X(v) = (1/v) {ln[(l+v)/(l-v)] - i7t}, s > 4ml. (40) 

The result (40) facilitates the sum-rule determination of closed loop contributions to Fq^. For 
example, the one-loop purely perturbative contribution [Fig. 4a] to the longitudinal component 
of the axial- vector current correlator (29) is given by [15] 

Cp[v] = (-3m^/27u2)[v2X(v) -i- divergent constant]. (41) 

The contribution of (41) to Fq j is easily obtained via a distortion of the contour C(Sq) to that in 
Fig. lb: 



[Fo(So)]p = (1/^)^ dsIm{CL[v]} 





s, 







= (3mJ/27C2)/ ds (1 - 4m2/s)'" 
4ml 



= (3m2/27C=^){so(l-4m2/so)^ + 2m^ ln|[l-(l-4mJ/so)^]/[l+(l-4m2/so)^] |} 

(42) 
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[F^(so)]p = (l/7c)/°dssIm{C^[v]} 





= (3m^/27u2) / ds s (1 - Am^sf 



4m^ 



= (3m2/47u2)(s2-2m2So)(l-4m2/so)^ 

+ (3m6/7i2) ln|[l-(l-4m2/sof]/[l+(l-4m2/so)^]|}. (43) 

We note that the results (42) and (43) are exact expressions obtained from the one-loop 
expression (41). To leading order in the quark mass m , one finds from (42) and (43) that 

[Fo(So)]p = (3m2so/27u2) + 0(mJ), (44) 
[F^(so)]p = Oml^yAn^) + 0{m\). (45) 



Two-Gluon Condensate Contributions to Fq j 



The OPE coefficient Cq2(Q2) is extracted from the OPE coefficient Eqi(Q^) in the 
"normal-ordered basis" [i.e., the "heavy quark" coefficients listed in Appendix B of ref. 15] as 

follows : 



Cg2(Q2) = Eg2(Q2) + (l/127l)C-q(Q2) - (mq/27l) In(m2/|i2) C^(Q2), (46) 



The linear combination (46) represents the coefficient in the "minimally-subtracted basis," which 
is chosen so as to avoid mass-singularities [16]. For example, in the normal-ordered basis 
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the coefficient of (cCsG^) for the longitudinal component of the axial-vector correlator is [15] 



E^2(Q^) = (-l/967uQ2)[16mj3+9v2)X(v)/(vV)] + [1/(487Uv''Q2)][9v^+ 4w^+ 3] 

(47) 

which generates an expansion 



E^2(Q^) = 1/37UQ2 - 5m^/67lQ4 + (mJ/7lQ^)[13/3 + 2 ln(mJ/Q2)] + 0(m^^Q^). 

(48) 

The leading term on the right hand side of (48) does not vanish as m^ — > 0, despite the chiral 
invariance of (afi^). Moreover, the right hand side has the quark mass appear in the logarithm, 
which could (in principle) lead to a large logarithm after subtractions. The change of basis (46) 
eliminates both problems, as is evident from direct substitution of (48), (34) and (30) into (46): 

Cg2(Q2) = - mJ/27uQ" + (mJ/7uQ6)[ll/3 - 2 hi(Q'/\i')] + 0(m^^Q\ (49) 

The result (49) is consistent with the general form (28), although the recipe (46) requires 
further modification if 0[mqln(mq/Q2 )/Q ] terms are to be eliminated. It is worth noting that the 
change of basis (46) differs from an operator redefinition proposed on chiral symmetry grounds 
in ref. [15] only by the presence of the final C]y[(Q2) term, which has already been shown not 
to affect the contour integrals leading to Fq and Fj. In Appendices B and C, the full contribution 
of C^(Q2) to the Fq 1 sum rules for the longitudinal component of the axial- vector correlator is 
determined to all orders in m^ by careful consideration of the C(sq) contour. However, 
contributions to Fq and Fj from the (a^G^) condensate can be evaluated to 0(mq) from 
application of (36) and the Cauchy residue theorem to (49): 

[F^(so)] = (l/27ti)<a,G ^ )/ ds C^2(-s) = (m^/Tcs^Xa.G ^ ) (50) 
<a,G^> t(so) 

[Fi(so)] = (l/27ui)<a,G2>[ ds s C^(-s) = [mj/27u + 2mJ/7USo]<a,G2>. (51) 

(ap^) C(So) 
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Higher-Dimensional Gluon Condensate Contributions to Fq j 

As in (46), the OPE coefficient Cq3(Q2) in the minimally-subtracted basis can be 
extracted from the coefficient Eq3(Q2) in the normal-ordered "heavy quark" basis [15,17]: 



Cg3(Q2) = Eg3(Q2) + [l/(3607im2)] C-^(Q^) + [l/(127imq)] Cm(Q^). (52) 



This change of basis once again eliminates leading-order mass-singularities. To see this, 

we demonstrate application of (52) to FESR's by once again considering the relevant 

contributions to the longitudinal component of the axial- vector current correlation function. The 

OPE coefficient E^3(Q2) is given by [15] 

E^3(Q^) = X(v)[7 + 23v2 + ISv'* + 5v^] 

247cQ^v^ 



+ [105 + 65v2 - 494v^ + 266w^ + 5v^ - 75w^% 

28 807cQV(1-v2) 



which generates the following expansion in inverse powers of Q ^ 



(53) 



E^3(Q2) = 



7tQ = 



a + 



90m^ 90Q2 45Q^ 



0(mj)] (54) 



The leading term of (54) diverges as m^^ 0, an explicit mass singularity. The next-to-leading 
term fails to vanish in the chiral limit. However, both of these terms, as well as the explicit 
0(mq) term in (54) cancel in (52) against corresponding terms from C-^ (30) and (34). 
Consequently, the OPE coefficient Cq3 is explicitly 0(mq), and is therefore suppressed relative 

to Cq2. 
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The suppression of Cq3 relative to Cqi in the operator product expansion appears to be 
a general property [17,18]. Corresponding factors of Cq2 and CqS for the scalar and vector 
current correlation functions, as well as for the transverse component of the axial-vector 
correlation function, as extracted via (46) and (52) from the "heavy quark" expressions in [15], 
are also seen to exhibit suppression by m^: 

By contrast, the dimension-8 contributions to scalar, pseudoscalar and vector correlation 
functions (which in our conventions are defined to have dimensions of mass squared) are shown 
in ref.[19] to be of the form [Aq + Bq ln(Q2/|x2)]<GVQ^, where \ and Bq are numerical: 
suppression by m^ does not seem to occur. For the longitudinal component of the axial-vector 
correlator, which picks up a factor of m^/Q ^ relative to the pseudoscalar correlator, a dimension- 
8 contribution to will then be proportional [via (36)] to mqBQ(G'^>/sQ. Such a contribution will 
be small compared to that of the dimension-4 condensate (G^> [eq. (51)] provided Bq(G^> is 
small compared to (G^>Sq, suggesting, in the absence of m^-suppression factors, that any further 
suppression of 2d-dimensional gluon condensates to FESR's is contingent upon the ratio (G^^)/Sq 
being small. Such a small ratio can be anticipated via dimensional and factorization arguments 
[e.g. (G^) < Sq]. 



V. Direct Single-Instanton Contributions to F\ 

In the instanton liquid model, the direct single-instanton contribution to the Rg Laplace 
sum rule (2) for the pseudoscalar (P) correlation function has been found to be [20] 

oo 

RP(t)= (1/7U)/ Im{[nP(s)]i„3j e-^Ms 





£P_ g-p V2t [K^(p , /2x) + Ki(p 2 /2x)], (55) 
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where p [= 1/(600 MeV)] is the instanton-size parameter. Since the pseudoscalar correlator is 
related to the longitudinal component (L) of the axial-vector correlator by 



U\s) = 4m2nP(s)/s, 



(56) 



L -4 



we see that the instanton contribution to the corresponding FESR Fj is 



^i(so) = (l/7t) / hn{[nL(s)]i„,J s ds 





(4m2/7t) / Im{[nP(s)]i„,J ds, 




(57) 



which is related via Laplace transformation to the function (55) for Rq(s) as follows: 



^/(t) =(4m2/7c) Im{[nP(t)]i„3j, 



m^i'm = I ^/(t) e-^' dt = 4m^ R^(s) = s " ^li^)- 





(58) 
(59) 



We see from (57) that ^^(0) = 0, and find that 



^i(t) = Sf-'[4m2 rP(s)/s] 



-1 



3p2m; 



c^-p^/2s [Kq(p2/2s) + Ki(pV2s)] 



(60) 



From (56), the direct single-instanton contributions to Fq and F^" are both 0(mj). This 
implies that the instanton contribution to Fq is small in comparison to (31), the leading 
(quark-condensate) contribution to Fq, which is why we are only concerned here with 
instanton contributions to FV. 
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Our use of the variables s and t is to retain consistency with standard Laplace transform 
conventions; the variable t will ultimately be identified with the continuum threshold Sq, and the 
variable s corresponds to the Borel parameter x in (55) [as defined in (2)]. The inverse transform 
of (60) may be obtained from the asymptotic expansions of Kq and Kj [21]: 

Ko(z) + Ki(z) = {nllzf e-^ [2 + l/(4z) - 3/(64z2) + 15/(512z3) ... ] (61) 



3m„^p 

m) = il 



4p2 32p4 



128p* 



(62) 



Using (62) and replacing t with Sq, we find that 



^i(So)=-^ G(2ps^), (63) 

where the function G(2pSo) is defined via 

G(w) = {[-w2/4 + 25/32 + 0(l/w2)] sin(w) + [-7w/8 + 15/(64w) + 0(l/w^)] cos(w)}. 

(64) 

The results (62-64) are not useful unless w 2pso > 1. Since Sq is generally expected to be 
at least 1 GeV, the expansion in large w is appropriate and useful [p"^ = 0.6 GeV]. In the large 
Sq limit, the leading perturbative contribution to F^" [eq.(45)] dominates the instanton contribution, 
which is at most linear in Sq (63-4). However, for values of Sq near 1 GeV 2, the instanton 
contribution is shown in the next section to be larger than the perturbative contribution, with 
phenomenological implications for the light quark mass. 

VI. Discussion: FESR's in the Pseudoscalar Channel and the Light-Quark Mass 



An old [22] and ongoing [23] controversy in sum rule applications concerns the 
failure of the field-theoretical content of the QCD sum rules to saturate the pseudoscalar channel. 

The essence of this problem is evident from a qualitative examination of the Rg and Laplace 
sum rules for the longitudinal component of the axial-vector current correlation function, as 
defined in (2) and (4). For suitable values of the Borel parameter x (M = » m^), one finds 
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that 



Ro= f^m^+ E F^M? expC-M^T) = -4<mqqq> + 0(m2), (65) 

M?< So 

a result consistent with the current-algebra GMOR relationship f^m^ = -4(mqqq> [24] as long 
as the subsequent subcontinuum resonances in the summation on the hadronic side of (65) are 
either sufficiently heavy (M? » 1/t » m^), or their decay constants F? are sufficiently small 
(F? « f^m^/Mj). The leading field- theoretical contribution to the Rj sum rule, however, is 
quadratic in the quark mass [1,10]: 

Ri= f^n4+ E F^M^ exp(-M?T) = m2[-4<mqqq> 3/(27u2t2) + <a^G2>/27U 

M?< So 

+ 4487ix<a,(qq) ^ >/27 + ...]. (66) 

Naively, the field-theoretical content of (66) is of order m^ times the field-theoretical content of 
(65), whereas the hadronic content of (66) is at least of order m^ times the hadronic content of 

(65) , suggesting that m^ and m^ are comparable. A thorough treatment of QCD contributions to 

(66) still yields substantially larger values of the light quark mass [10-12] than are anticipated 
from other phenomenology [5], as already noted in the Introduction to this paper. 

This mismatch in scale [i.e., Rj/Rq ~ m^; R^^^/Rg*"^^ ~ m^] superficially characterizes the 
FESR's Fq and F^ as well. However, these FESR's provide a much cleaner framework for 
extracting limits on m^, enabling one to avoid the large-width modifications to the hadronic- 
resonance content of (66), as discussed in Sections 11 and HI, as well as higher-dimensional 
condensate contributions (including that of (a/^qq)^}) to the field-theoretical content of (66), as 
is also discussed in Sections IV and V. 

Direct single-instanton contributions (55) to the Laplace sum rule have been argued in a 
number of places [1,20,23] to be necessary for the saturation of the pseudoscalar channel. If we 
incorporate such contributions (63,64) into the FESR F^, in conjunction with the (width- 
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independent) hadronic contributions (19) as well as the leading [0(mq)] field theoretical 
contributions (32,45,51), we find that 



F\= fX+ ^ FlMt = m^[-4(mqqq> + (a,G2>/27U + 3s^47u2 

M?< Sn 



+ (3/ny) G(2ps^) + 0(mq)]. (67) 



For each subcontinuum pion-excitation state, we define the parameter 

r, = (FjMt)/(flm'^. (68) 

We can then rearrange (67) to obtain the following relationship for the Ught-quark mass: 

ml = _imliL+_^_£i)_ , (69) 
A + [G(w) + w^/64]B 

where... 

1) ... the summation is understood to be over only those resonance peaks below the 
continuum threshold (M? < Sg); 

2) ... the dependence on the continuum- threshold Sq enters through the variable 

w = 2ps^; (70) 

3) ... the function G(w) is given by (64); and 

4) ... the constants A and B are given by 

A = -4<mqqq> + <a3GO/27C, (71) 

B = 3/(%Y)- (72) 

The relationship (69) should retain approximate validity as Sq increases to include (in full, as 
noted in Section III) the resonance peaks of additional pion-excitation states. In particular, one 
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would expect the contribution from 11(1300), the first pion-excitation (M = 1300 ±100 MeV, F 
= 200 - 600 MeV [5]) to be fully subcontinuum if Sq > 4 GeV^. Possible additional 
contributions may accrue in full from 11(1770) and X(1830) at even larger values of Sq. 

Using standard parameter values [(mq'qq) = -f^m^4, (afi^) = 0.045 GeV"^, p"^ = 600 
MeV], one can then estimate the following numerical lower bound on the quark mass from (69): 



m = |x(w) Vl + E r^ > |x(w) Vl + r^@(sQ- 4 GeV^), (73) 



, , 2.6 MeV 
|i.(w) = 

{0.0075 + 0.039[G(w) + w4/64]}'' ' ^'^^ 



where r^ is just the value of (68) appropriate for the first pion-excitation state. Although chiral 
Lagrangian arguments have been recently advanced suggesting that rj is substantially less than 
unity [4], sum-rule estimates for rj of order unity and larger [13] have received further support 
[25] from recent Laplace sum-rule fits. 

We reiterate that the FESR-based inequality (73) [and the relation (69) from which it is 
derived] avoids any need for a narrow-resonance approximation, which would certainly be 
unphysical for dealing with broad subcontinuum pion-resonance states. The QCD-vacuum 
condensates that contribute are all lumped into the constant A (71); condensates such as 
(qG-aq), {a^(qq)^), and (a^G'') do not generate any 0(m^ contributions to as has already 
been discussed in Sections III and IV. Even dimension-8 gluonic condensate contributions can 
be expected to be suppressed relative to those of (a^G^) by the dimensional arguments presented 
at the end of Section IV. 

In Table I, we tabulate the quark-mass lower bound |i(w) for values of Sg ranging from 
1 GeV 2 to 4 GeV ^ . We also tabulate the same function in the absence of instanton 
contributions [i.e., with G(w) = 0] in order to demonstrate the key role instantons play in 
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obtaining a lighter and phenomenologically consistent quark mass over the entire range of Sg 
considered. When the contribution of instantons is absent (Column 4 of Table I), we find that 
|x(w) decreases from 9.1 MeV by a factor of four as Sq increases from 1 GeV^ to 4 GeV^ . This 
behaviour, if taken seriously, would not only suggest via (73) a rather large quark mass (~ 9 
MeV), but also a very large aggregate contribution 

Erj » 15 (75) 

from subcontinuum resonance-peaks as Sg increases to 4 GeV ^ . 

The instanton term G(w) in the denominator of (74) greatly ameliorates these effects. 
When the instanton term is included (Column 3 of Table 1), we find that |l(w) decreases from 
5.7 MeV by only a factor of two as Sq goes from 1 GeV^ to 4 GeV 2, suggesting via (73) a 
lighter (~ 6 MeV) quark mass in conjunction with a phenomenologically reasonable aggregate 
contribution 



(76) 



from subcontinuum resonance-peaks as Sq increases to 4 GeV ^ . In view of the sparseness of 
such pion-resonance states [which suggests replacing E rj with rj], it is noteworthy that the 
estimate (76) is quite compatible with past and present sum rule estimates for rj [13,25]. 

It is best to regard the results presented in this section as essentially qualitative. We have 
utilized only the one -loop-order purely-perturbative contribution to the correlation function 11^ 
- higher-order terms can be expected to alter the coefficient of the w'^-dependence in the 
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denominator of (74).^ The key point here, however, is that the function G(w) arising from 
instantons is oscillatory (64), going from positive to negative values as Sq increases from 1 GeV ^ 
to 4 GeV 2 . Moreover, G(w) is not only positive, but is also larger over the range 1 GeV ^ < 
Sq < 1.6 GeV 2 than the factor w'^/64 arising from perturbation theory, thereby lowering and 
stabilizing the quark mass (73) in a region for which there is at most only a partial contribution 
from the lowest subcontinuum resonance. 

Acknowledgments : 

ASD, VE and TGS are grateful for support from the Natural Sciences and Engineering 
Research Council of Canada. VE is also grateful for several useful discussions with V. A. 
Miransky. AHF and YX wish to acknowledge research funding from the Faculty of Graduate 
Studies of the University of Western Ontario. 



Inclusion of the renormalization group (RG) dependence of the running quark mass also 
lowers somewhat the size of the aggregate resonance contribution Er^ as Sg — > 4 GeV 2. 
Assuming Aq(-.q = 0.2 GeV, we find near-constancy of the RG-invariant quark mass m {m^ — > 
mq(so) = m/[ln(Vs'o/AQCD)]^^'^} over the 1 GeV 2 < Sq < 4 GeV^ range of Table I provided Er^ 
2 if instantons are included, with Erj — > 10 if instantons are not included. 
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Appendix A: Imaginary Parts of Correlation Functions for s > 4m^ 



The purely perturbative contribution to the longitudinal component 11^ of the axial vector 
correlation function, defined via 

ijd^x eiP''<0 Tj^5(x)j^g(0) 0> 

= [g^iv-PuPv/p'^n^Cp') + [p^pyp2]nL(p^) (A.1) 



with j|^5(x) = u(x)Y|^Y5 d(x) can be obtained from Fig 4a: 



[nL(p')]pen-Cp,>') 



3m- 



n-4 



Ye + In 



i(p') 



(A.2) 



J 

i(P^)= / 



dx In 



1 - i_x(l-x) - i 8 



(A.3) 



m' 



If p ^ < 4m ^ , the argument of the logarithm is positive definite, and the i 8 factor is irrelevant 
to the evaluation of the integral. It is straightforward to find that Cp^^^ is real provided 

p^ < 4m ^, and that 
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I(p2) = -2 + 



4m 2 



In 




N 



4m ^ 



, < 



(A.4) 



I(p2) = -2 + 2 



4m 2 



1 tan 1 



N P' 



4m 2 



Vl/2 



< p < 4m' 



(A.5) 



One easily finds from either expression that Lim I(p ) = , and from the latter expression that 
Lim I(p ^) = -2 . The results (A.2) and (A.4) are consistent with C ^ as calculated in ref. 

p2^(4m2)- 

[15]. Utilizing the notation of that reference, one sees that 



X(v) = 1 In 

V 



v + 1 



v-1 



_ I(p')+2 



(A.6) 



with 



V = 



4m 2 



(A.7) 



The relationship (A.6) between X(v) and I(p ^), the latter quantity defined via the integral (A.3), 

can be utilized to determine the imaginary part of X(v) when p^ > 4m ^, as shown below. 

If p^ > 4m ^, the argument of the logarithm in the integrand of (A.3) can be factorized 
and integrated by parts as follows: 
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with 



and with 



I(p') = In 



J dx ln(x - - ie^' ) 



1 

J dx ln(x - X + ie'^ ) 



= -2 - (X, + ie^O J 



dx 



X - X^ - 18 



(X - ie'O J 



dx 



X - X + 18 



(A.8) 



1 ± 



x+ = 



4m 2 



N P' 



(A.9) 



,// _ P I e 



m\x^ - x_) 



(A.10) 



Note that if > 4m ^, then 0<x<x^<l, and, hence, that e" > 0. Consequently, the 
pole in (A.8) at x^ + ie^' is above the real x axis, and the pole at x - ie'^ is below the real 
X axis, permitting the equivalent contours of Fig. 5 to run below x^ and above x . Using the 
contours of Fig 5, with C+ and C assumed to be semicircles of radius 5 about x = x^ and x , 

10 

respectively, one finds that [C+: z = x+ + 5e -; range of 0+: TC — > 2%; range of : 7C — > 0] 
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1 

/ 



T -5 



dx 



- Lim 



+ + 



dz 



In + i7t , 



(A.11) 



1 

/ 



dx 



X-T +18 



In ^ - i7U 



(A.12) 



Substituting (A. 11) and (A.12) into (A. 8), one finds for > 4m^ that 



I(p^) = -2 + w 



In 



^ 1+v ^ 



17t 



(A.13) 



Using the relationship X = (I+2)/v , we then see that X develops a negative imaginary 



part: 



ImX ^ -— , p^>4m- 

V 



(A.14) 



Appendix B: Evaluation of the Gluon Condensate Contribution to Fq 

The "heavy-quark" (h.q.) gluon condensate contribution to Ilj^, as defined in (A.l), is 
obtained from Appendix B.3 of ref [15] as the sum of coefficients [Ciq2]^ and [C2q2]^ for 



the axial-vector current correlation function 



s = p^, V = Vl-4m^/s 
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a 



[3(l-v2)2x(v)-6(Uv2)] , 



a 



1- 967tsv^ 



[3(l-v^)2 (1+v^) X(v)-2(3-2v^+3v'')] 



[CiG^ - C^G.]^ - ccE^. = aEp^j^ + aC^ X(v) 



oE, 



pole 



a 

967t 



18 14 

+ 



24m- 



2^2 



s-4m^ (s-4m'^) 



111 



1 3 

+ 



3 4 3 2 

S V S V 



We have extracted a factor of a so that Eq2 as defined in (B.4) is consistent with Eq 
in Sec. 3. 

The gluon condensate contribution to the finite energy sum rules 



Fo^ = ^ / n,(s) 



ds 



C(s„) 
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Fi"^ = J Hl (s) s ds , (B.8) 

C(So) 



can be obtained via (46) and (50) from direct evaluation of the integrals 

Gq= J Eqz ds , (B.9) 

C(So) 

Gi = J Eq2 s ds , (B.lO) 

C(So) 

with the contour C(Sq) distorted as in Fig. 6 to encompass any pole singularities of Eqq at 
s = or 4m ^ as well as the branch singularity for s > 4m ^. Using Eq. (A. 14), one finds that 



lin f ^ 

. < V 



Gg = -2i% I _ ds 



4m 

+ J Epoie ds + J Ep^i^ ds 

(B.11) 



+ 



J X(v) ds + J X(v) ds , 



where the contours Cq and C^^^i are clockwise circles of radius E about s = and s = 4m ^, 
respectively. We see from (B.5) that 
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c 8 



J Ep„ie ds = ^ . (B.13) 



4m 



The remaining three integrals in (B.ll) are evaluated as follows. Using the expression 
for in (B.6), we find that 

Sq 

-2i% J ^ ds = -i m^ [I^ + 3I2] , 

4m 

where 

Ii = ^ ds = f s-'%-4m'r''^ ds , (B-15) 

4m 4m 



r°^ds= [' s-3/2(^-4mVds . (B-16) 

4m 4m 



Both integrals can be evaluated via the trigonometric substitution s = 4m ^sec^S . One then finds 
that 



8m 4 i 



cos 
sin^e 



d0 



1 



24m 4 
1 



sin^e 



0.. 



8m' 



sin0 



0. 



(B.17) 
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I. = 



8m ^ i 



cos 6 
sin^e 



de 



8m' 



sine 



+ sine 



(B.18) 



where, using the parameterization of (A.7), we find that 



sec 



-1 



1/2 



sin e.. = \ I -4m^/s, 



2/„ _ 



(B.19) 



'0 



and that 



sec 



(4m ^+8) 



1/2 



2m 



sine. 



4m ^+8 



1/2 



(B.20) 



Substituting (B.19) and (B.20) into (B.17) and (B.18), we find from (B.14) that 



-2i7l 



J V 



ds 



4m 



1 

"8 



3Vn 



2 o 

- 3Vn 



im Sim 



+ 0(8^^2) 



3e^'^ 88 



1/2 



(B.21) 



The integral around the origin is straightforward to obtain from (B.6) and (A. 6). The 
integrand 
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C. X(v) = 



1671 



2+I(s) 



-3/2 



6m 



8m ^ 



s-4m2 (s -4m (s -4m 2)3 2s 



(B.22) 



has a simple pole at s = because 1(0) = [as noted below (A.5)]: 



J X(v) 



ds = - — 
8 



(B.23) 



Note that (B.23) exactly cancels (B.12), indicating that the origin can be excised from the contour 
of Fig. 6. 



This cancellation is not peculiar to the channel we are in. We have verified (Appendix 
D) that an identical cancellation occurs in the scalar, vector, and transverse-axial channels 
between the contributions of explicit s = poles in Eq2 [as in (B.12)] and the integrals of 

Cjj^X(v) portions of Eq2 around Cq [as in (B.23)]. Thus the quantum-field-theoretical 

singularities in Gq and Gj all occur for s > 4m 2 on the real s-axis for all of the above- 
mentioned channels. 

The divergence as 8 — > in (B.21) is cancelled exactly by the integration of C^^^ X(v), 

as a given in (B.22) over the contour C^^i around s = 4m 2, a cancellation which also occurs 
in the other three channels mentioned above. This cancellation is most easily seen by continuing 
the expression (A.5) to complex values of s in the vicinity of s = 4m 2 : 
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( ^l/2 



I(s) +2=2 [(4m2-s)/s]^^^ tan ^ 



7U 





1/2 








4m ^-s 


- 2 


4m ^-s 


2 

+ 


4m ^-s 


I ^ J 




I ^ J 


3 


I ^ J 



(B.24) 



On the contour s = 4m ^ + 8e^® with a clockwise rotation of 9 from 27U to 0. When 

s > 4m ^, the correct (negative) sign of the imaginary part 2i Im {l(s)} = I(s +i 5 ) - I(s -i 5 ) 
is obtained by requiring that 

(4m2-s)i/2 ^ _-gi/2 gi9/2 ^ (B.25) 



as 



2i Im {l(s)} = Lim 

9^0 



71 



f ^1/2 

4m ^-s 



V J 



Lim 
e->23t 



71 



( Al/2 

4m ^-s 



V J 



(B.26) 



with s = s(e) = 4m^ + 8 e'®. Upon substitution of (B.24) into (B.22) one finds that 



J X(v) 



ds^ 



7i 
24 



16 



. J s -1^^ (4m2-s)-i^2 

C. 2 
4m 

6m 2 J s (4jjj2_g) 3/2 
8m ^ J s (4m2-s) ^^2 ds 



^4012 



J (4m2-s)^^2 ds 



C4m2 



(B.27) 
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The factor -7i/24 is just -2%i times the aggregate residue at s = 4m ^ obtained from 

multiplication of (B.24)'s integer powers of (4m ^-s) into (B.22). This pole contribution 
explicitly cancels the pole contribution (B.13). The remaining integrals in (B.27) result from 

multiplying the leading 7c[(4m^ -s)/s]^^^ term of (B.24) into (B.22). These integrals are easily 

evaluated around the clockwise contour C^^i via (B.25): 



J (4m^-sy^'^ s '^'^ ds = 0(e^'^) 



4m 



r (4m 2 -s) -3/2 s ^ 2l ^ ^ ^^^^^ 



(B.28) 



f (4m2-s)-5/2 s ^ ^ 1 ^ o(8i/2) , (B.30) 



J (4m2-s)1^2 g 3/2 ^ Q(g3/2) 



c , 24 881^2 



(B.31) 



Substituting (B.28 - B.31) into (B.27) we find that 

f C, X(v) ds ^ + ^ + + 0(8i^2) (g 32) 



explicitly cancelling the divergencies in (B.21). Since all the s=0 and s=4m^ pole terms 
contributing to Gq have also been shown to cancel, we find that Gq is equal to the upper-bound 
contribution of the first integral on the right-hand side of (B.ll): 
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Go = 



-2m 



J V 



ds = ± 
8 



+ — + 3v. 



(B.33) 



To obtain the full contribution of < a^G ^ > to the Fq sum rule, we substitute Eq. (46) from the 
text into (50), utilizing the results (B.33) in conjunction with Eqs. (30) and (34) from the text: 



J <a.G^> s 



16k 



m 



1 2 , 

— + + 3Vn 



3v.^ 



14m' 



1 

37t 



TUSn 37USn 



(B.34) 



Appendix C: The Gluon Condensate Contribution to Fj 

Consider first the integral Gj (B.IO), which can be evaluated via the following integrals 
arising from the distortion of C(Sq) indicated in Fig 6: 

Gj = -2iK J s ds + J Ep^j^ s ds 



jEpoieSds . JC, X(v) 

--4m 2 ''O 

J C^ X(v) s ds . 



s ds 



(CI) 



One sees from (B.5) that 
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J Epole s ds - , 
Co 

J Epole S ds = . 

Using the expression for C^^^ in (B.6), we find that 

C 

-2in J ^ s ds = -im'^ [I3+3I4 

4m^+e 



where the integrals I3 and are evaluated using (B. 19-20), as in the previous section 



"0 



4m ^^+6 



- ds = - J_ 

5 6m^ 



sin-^0 



4m 1 i/2s 



"0 

/ 4- 

_2.„ S^V 



4m-' +E 



3 2m2 



e.. 



sinG 



2m \q me^/2 



Substituting (C.5) and (C.6) into (C.4) we find that 
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Using (B.22), we find that X(v) s has no poles at s = [note that 2 +1(0) = 2], in 
which case 



J X(v) 



s ds = 



(C.8) 



Once again, we note that the origin can be excised entirely from the contour of Fig 6. We have 
verified explicitly that integrals (C.2) and (C.8) are zero in the scalar, vector and transverse axial 
channels as well (see Appendix D). 

As in the previous section, the divergence in (C.7) as 8 —> is exactly cancelled by 
integration of C^^^ X(v) s around the contour C^^i. From (B.22) we find that 



C, X(v) s = 



- [2+I(s) 



m 



m 



(s- 4m 2)2 (s-4m^) 



2x3 



(C.9) 



If we substitute (B.24) into (C.9) and integrate around C^^2, we easily separate a pure-pole 
contribution from an 8 -dependent contribution involving half-integral powers of (4m ^-s): 

J C^ X(v) s ds = + 2m* J (4m^-s)^'^ s '^'^ ds 

C4m2 ^4m^ 

- 2m^ J (4m2-s)"^'2 5 -1/2 



_ 5im2 ^ 7im3 ^ 4im^ ^ . (C.IO) 

3 2e"^ 383^2 

The final line of (C.IO) is obtained through use of (B.29) and (B.30). Not only are the 8- 
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dependent terms in (C.7) cancelled by the final line of (C.IO), but the pure-pole contribution 
(C.3) also cancels against the pole term in (C.IO). Thus we find, as in Appendix B, that Gj is 
equal to the upper-bound contribution of the first integral on the right-hand side of (C.l): 



Gj = -2i7t 



r c s 

ds 

J V 



im' 



1 



6Vn 



2vp 



(C.11) 



To obtain the full contribution of <cc^G ^> to the Fj sum rule, we again substitute Eq. (46) of 
the text into (51), utilizing the results (C. 1 1) in conjunction with Eqs. (30) and (34) from the text: 



Fi^so) 



1 3 

+ 



6Vn 



2Vn 



2 
3" 



m 



<a G^> ■ 



1 + 



4m' 



14m* 160m 6 

+ 

2 ™ 3 

Sn jSn 



(C.12) 



Appendix D: Gluon Condensate Contributions to Fq j in Other Channels 

Utilizing the notation and conventions of Appendices B and C, we list the following 
results obtained from scalar, vector and the transverse component of the axial-vector correlation 
functions: 

Scalar Channel 
From Appendix B.l of ref 15, we have 
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We find that 



fo^U. " = •^(^Pole " X(V); 



(D.l) 



E - (3-v') 

pole 2 ' 

167tSV 



(D.2) 



^ _ (l-v^)(3+v^) 
327USV ^ 



(D.3) 



Epole ds = -i 



(D.4) 



(D.5) 



-2i7t 



< V 



4m^ + e 



ds = -L 
8 



3 

Vn 



6m 

7?. 



(D.6) 



J X(v) ds 



(D.7) 



J X(v) ds 



3i 



1 _ 2m 



4m 



(D.8) 



Summing (D.3-8) we obtain 
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/ 

C(s„) 



Eq2 ds 



1 

"8 



+ Vn 



(D.9) 



We also find from Appendix B.l of ref. 15 that 



J ds = 6i7U , 

C(So) 



(D.IO) 



J ds = , 

C(So) 



which implies via (46) and (60) that 



167t 



(D.ll) 



(D.12) 



Unlike the case of Fq, the FESR Fj requires the use of (46) to eliminate a logarithmic mass 
singularity in Gj , obtained by summing the following five integrals: 



/ EpoieS ds = , 



(D.13) 



J 3im ^ 
E„„i„s ds = 

pole 2 



(D.14) 



-2i7U 



^ V 



s ds 



im 



4m + E 



3 /i 2s 6m 
— + 2fn(l -Vq) + 

^0 



(D.15) 
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J X(v) s ds = , 



(D.16) 



J X(v) s ds 



Sim 3m " 



(D.17) 



We then find that 



C(So) 



Eq2 s ds = 



im 



-_ + 2Cn 



4m 2 



(D.18) 



which is not analytic in m at m = . However the results 



J C-q s ds = 4im27U , 

C(So) 



(D.19) 



J C]y[ s ds = -2im7U 

C(So) 



(D.20) 



used in conjunction with (46) and (51) eliminates the quark-mass from the logarithm: 



Fi(so) 



m 



<«sG'> 2% 



( \ 



3 1 

+ 

3 



2v, 



to 







4)1 



<a^G2> 



(D.21) 



Transverse Axial Channel 



From Appendix B.3 of ref. 15, we have 
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We then find that 



[ClGi.q. - « = «(Epole - X(V)) , 



(D.22) 



^pole 



Stusv^ 



(D.23) 



2^2 



167tSV^ 



(D.24) 



Epole ds = - 



(D.25) 



J 



Epole ds = - 



(D.26) 



-2i7l 



•0 ^ 

J - 

^' V 



ds 



4m + e 



1 2m 



(D.27) 



J C, X(v) ds 



(D.28) 



X(v) ds = - + 
4 



(D.29) 
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As before, the contour-radius singularity as 8 ^ cancels between (D.27) and (D.29): 



Go = J ds = 

C(s„) 



(D.30) 



Since in this channel, one finds that [15] 



J C-qq ds = -47Ui , 

C(So) 



(D.31) 



J Cm ds = , 

C(So) 



(D.32) 



we find via (46) and (60) that 



Fo(So) 



1 4 
Vq 3 



(D.33) 



Corresponding results for are listed below: 



Epoie s ds = 



(D.34) 



J Epoie s ds = im^ , 



(D.35) 
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-2m 



r — i s ds = im^ 
J V 



4m + 8 



1 2m 



J X(v) s ds = , 



/ 



X(v) s ds = -im^ + 



J 

C(s„) 



„ , 8i7tm ^ 



/ 

C(So) 



Cm s ds = 



m m" 



27tVn 9n 



Vector Channel 
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From eq. (11.19) of ref. 15, we find that 



_ _(3-2v2+3v4) 
487tsv^ 



''pole 



(D.42) 



^ ^ (i-vy(i+v^) 

327USV* 



(D.43) 



We then find Fq (D.51) from Gq, the sum of (D.44-48), in conjunction with (46) and (D.49-50): 



Epole ds = - 



(D.44) 



J 

4m 



(D.45) 



-2i7U 



s' V 



4m + e 



ds = _ 
8 



^0 + 



1 8m m 



3v' 383^2 gl/2 



(D.46) 



J 



Cx X(v) ds = -A , 



(D.47) 



J X(v) ds 



1 im im 

+ + 



4m 



24 3e3/2 ggi/2 



(D.48) 



J C-qq ds = -4i7i , 

C(So) 



(D.49) 
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J Cm ds = , 

C(So) 



(D.50) 



[Fo(so)] 



1 



^0 + 



3Vn 



(D.51) 



Corresponding results for are listed below: 



/ Epoie s ds = 



(D.52) 



J 



c J 2im^ 
E ,^ s ds = 

pole 2 



(D.53) 



-2i7U 



'0 p, 

J - 

:/ V 



s ds 



im 



4m-^ + B 



1 1 

+ 



8m 3m 

3/2 el/2 



3v' 38^'^ £ 



(D.54) 



J X(v) s ds = , 



(D.55) 
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/ 



„ , J 2im^ 4im^ 3im^ 

C X(v) s ds = - + + 

3 3e3/2 2^"^ 



(D.56) 



/ 

C(So) 



J 16-2 

C- s ds = - — iTum , 
qq 3 



(D.57) 



J 

C(s„) 



Cm s ds = 



(D.58) 



m' 



47t 



1 1 
+ 



Vq 3v, 



8 

3 "9 





(D.59) 
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Figure Captions: 



Figure 1: a) The contour C(Sq). 

b) Distortion of C(Sq) to enclose the positive real s-axis. 

Figure 2: The Breit-Wigner resonance shape y(s) = mr/[(s-m2)2+m2r2] expressed as a 
sum of symmetric square pulses. For an infinitesimally thin pulse at a given value 
of y, the pulse will extend from s = m^ - [mP/y - m2r2]^tos = m2 + [mP/y - 
m2p2]^. If there are n such pulses, the j* pulse is at y = j/(nmP). 

Figure 3: a) Leading (qq) contribution to current correlation functions. 

b) Typical leading (qG-oq) contribution to current correlation functions. 

c) Typical leading {a^(qq)^) contribution to current correlation functions. 

Figure 4: a) Leading purely-perturbative contribution to current correlation functions. 

b) Typical leading (a^G ^ > contributions to current correlation functions. 

c) Typical leading (G'^) contributions to current correlation functions. 

Figure 5: Distortion of the integration contour along the real s-axis consistent with the 
location of the T+ singularities in the complex s-plane. 

Figure 6: Distortion of the C(sq) contour [Fig. (la)] for (a^G^ contributions to Fqj sum 
rules. 
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w 



So in GeV2 ^i(w) in MeV |i(w) in MeV 

[Eq.(70)] [Eq.(74), with G(w) given by (64)] [Eq.(74), with G(w)=0] 



3.3 


0.98 


5.7 


9.1 


3.5 


1.10 


5.2 


8.2 


3.7 


1.23 


4.8 


7.4 


3.9 


1.37 


4.5 


6.7 


4.1 


1.51 


4.2 


6.1 


4.3 


1.66 


4.0 


5.5 


4.5 


1.82 


3.8 


5.1 


4.7 


1.99 


3.7 


4.7 


4.9 


2.16 


3.5 


4.3 


5.1 


2.34 


3.4 


4.0 


5.3 


2.53 


3.4 


3.7 


5.5 


2.72 


3.3 


3.4 


5.7 


2.92 


3.2 


3.2 


5.9 


3.13 


3.1 


3.0 


6.1 


3.34 


3.05 


2.8 


6.3 


3.57 


3.0 


2.6 


6.5 


3.80 


2.9 


2.5 


6.7 


4.04 


2.8 


2.3 



Table I: Behaviour of |x(w) with increasing Sg in the presence (Column 3) and in the absence 
(Column 4) of direct single instanton contributions to the finite energy sum rule. 
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